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Abstract. We consider the class of Beltrami fields (eigenfields of the 
curl operator) on three-dimensional Riemannian solid tori: such vector 
fields arise as steady incompressible inviscid fluids and plasmas. Using 
^0 ' techniques from contact geometry, we construct an integer-valued index 

, for detecting closed orbits in the flow which are topologically inessential 

(they have winding number zero with respect to the solid torus). This 
| index is independent of the Riemannian structure, and is computable 

d • entirely from a C 1 approximation to the vector field on any meridional 

disc of the solid torus. 
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1. Introduction and summary 



Consider the class of Beltrami fields — the volume-preserving eigenfields of 
the curl operator. Such vector fields are the source of numerous interest- 
ing phenomena in inviscid fluids and plasmas. For example, Beltrami fields 
are the only steady three-dimensional Euler flows which admit chaotic La- 
grangian dynamics. Beltrami fields are also common approximations to the 
magnetic field lines in large-scale structures within the solar corona. 



a 

Despite their importance and inherent intricacy, very little is known about 
■ the dynamics of Beltrami fields apart from numerical simulation 22, F7| and 

Melnikov analyses of near-integrable Beltrami fields |18, |35], ^l], §T26| - an 
important but extremely small class of solutions. We consider the subtle 
problem of understanding how much and what kinds of dynamics Beltrami 
fields are forced to possess given the underlying topological features of the 
fluid domain. 



In a series of papers [13, [14], [15]], the authors develop techniques for deter- 
mining forced behaviors in steady inviscid fluids via the topology of contact 
structures, the odd-dimensional analogues of symplectic structures (see, e.g., 



[|], 11] for an introduction to contact geometry). In this paper, we give an 
application of these techniques to Beltrami fields on solid tori. One of the 
features of our topological approach is that it is independent of the Rie- 
mannian metric and is furthermore robust with respect to perturbations of 
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the vector field, without resorting to any hyperbolicity or nondegeneracy 
assumptions usually required to preserve closed orbits. 

We restrict attention to Beltrami fields on Riemannian solid tori, such as 
would occur in the case of a force-free plasma in a containment device. 



In [14], it is shown using techniques from contact topology and pseudo- 
holomorphic curves that steady Euler fields on an invariant solid torus always 
possess a closed orbit, independent of the Riemannian structure and volume 
form: 

Theorem: ( ||l4|| ) Any steady real-analytic solution to the Euler equations 



(2.i) on any invariant Riemannian solid torus possesses a closed orbit. 



Since every Beltrami field is a steady Euler field, the result holds true for 
all Beltrami fields. The restrictive smoothness assumption is necessary for 
using singularity theory arguments for the integrable Euler fields — in the 
setting of pure Beltrami fields, the techniques are valid up to smoothness 
class C 2 Wc 



In this paper, we define an integer-valued index for detecting the presence 
of contractible closed orbits — those closed orbits which can be shrunk to 
a point within the solid torus. Thus, we are not primarily concerned with 
the class of Beltrami fields which possess a cross-section (and hence trivially 
have a closed non-contractible orbit by the Brouwer fixed point theorem). In 
cases where no section exists, it is very difficult to determine the existence of 
closed orbits. Indeed, the recent examples of fixed-point-free vector fields on 
a solid torus without any periodic orbits (constructible via [^] in the real- 
analytic case, and by |28| in the C 1 volume-preserving case), demonstrate 
the delicacy of the problem. 

A related scenario to which our results apply is that of a Beltrami field on 
a long tube D 2 x IR which is periodic in the third variable. The problem of 
finding contractible closed orbits in the solid torus obtained by quotienting 
out the periodicity is precisely the problem of finding an orbit which is closed 
in the long tube. 

The index we construct is a type of linking number with respect to a con- 
tact structure — the so-called self-linking number of a transverse knot, well- 
known to contact topologists. The data required to compute the index is 
minimal: one needs information about the Beltrami field along some (arbi- 
trary) meridional disc in the solid torus. The principal contribution of this 
note is to retool the contact-topological index in a form which requires no 
knowledge of the contact structure per se. Since contact structures are very 
stable with respect to perturbation, the vector field need only be known 
approximately (C 1 ) along the disc. Hence, this index can be computed 
numerically with full rigor. 
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Sections ^ through || assemble the relevant ingredients. Section |5| presents 
the technical result on contact structures used to define the index and prove 
its major properties. The section §|| gives a very simple method for comput- 
ing this index from a minimal amount of data: one need simply know what 
the vector field X approximately looks like near some finite number of points 
on a meridional disc of the solid torus. It is our hope that this computability 
may allow for utilization of this index in the analysis of experimental data. 



2. Beltrami fields on Riemannian manifolds 



Let M be an arbitrary 3-manifold with Riemannian metric g and (arbitrary) 
volume form p,. Given a vector field X on M, one can consider the dual 
I-form g(X, •) to X that pairs with a vector Y to give the inner product 
g(X, Y). In this general setting, the curl of a vector field X on M is the 
unique vector field V x X satisfying 

(2.1) xX,-,-)=d(g(X,-)), 

where d denotes the exterior derivative on forms. The curl operator is linear 
and its /^-preserving eigenfields are known as the Beltrami fields. In other 
words, X is Beltrami if and only if it is volume-preserving and V x X = XX 
for some constant A. One can also consider the class of eigenfields with 
scalar fields as eigenvalues: V x X = fX for / : M — > K. Our techniques 
are adaptable to this more general format, but we restrict to pure eigenfields 
here for simplicity. 

Beltrami fields arise in several contexts: 



1. Beltrami fields are always steady solutions to the Euler equations for 
an inviscid incompressible fluid 

du 

(2.2) — + V u u=-Vp ; C u u = 0, 

where V u u is the covariant derivative of the velocity field u along itself, 

1 1 1 2 

and p : M — > R, the pressure function, can be chosen to be ^ \\ u \\ ■ The 
Lie derivative C u fx of the volume form along u vanishing is equivalent 
to u being divergence-free. 

2. Beltrami fields also yield steady solutions to the ideal MHD equations 

du 

— + V u u = - Vp + (V x B) x B 

(2.3) ^-Vx(ux8) = 

C u fi = CbH = 0, 

where B denotes the magnetic field. In this context, Beltrami fields 
are known as force-free fields. 
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3. Beltrami fields are all extrema of the 1? energy functional 



(2.4) ||u|| 2 := £ / \\u\\ 2 dn 

1 JM 

under the action of the volume-preserving diffeomorphism group of 
M. Eigenfields of curl having the smallest nonzero eigenvalue globally 
minimize the energy |2|, ||. 



Beltrami fields also play a role in the analysis of the stability of matter 1 3C ] 
and in the formation of dynamos ||. 

The topology and dynamics of Beltrami fields are subtle: witness the com- 
plex dynamics of the well-known ABC fields on the Euclidean 3-torus 
The existence of fixed-point-free Beltrami fields on general Riemannian 3- 
manifolds is highly nontrivial [O] , as is the presence of closed orbits within 



such fields PL 14 



3. Contact structures and topology 



Contact structures are the natural complements to Beltrami fields. Loosely 
put, a contact structure on a odd-dimensional manifold M is a hyper- 
plane field which is maximally nonintegrable. More specifically, on a three- 
dimensional manifold, a contact structure is a smoothly- varying plane field 
(a choice of a two-dimensional subspace £ p in each tangent space T p M) which 
cannot be stitched together into leaves of a foliation, not even at a point. 
Locally, every contact structure £ is the kernel of a differential 1-form a 
satisfying the contact condition: 

(3.1) a A da ^ 0. 

Otherwise said, a Ada is locally a volume form on M. Any 1-form satisfying 
( |3.1|) is called a contact form. If a A da is a globally defined volume form on 
M, then the contact structure is said to be cooriented: all contact structures 
which arise in connection with Beltrami fields are of necessity cooriented, 
and we will restrict entirely to this category. 

A contact form a on an oriented three-manifold M is said to be positive if 
the sign of a Ada is positive with respect to the orientation of M . Otherwise, 
a is said to be negative. The sign is a property of the contact structure and 
is independent of the defining 1-form. 

Canonical examples of positive contact forms on M 3 include dz + x dy, dz + 
r 2 sinrd6, and cosrdz + r sin rd6, the latter two being given in cylindrical 
coordinates. 

Much of the current interest in contact structures arises from fairly recent 
elucidations of their topological and dynamical properties (see, e.g., jj], |i~i"||). 
Studying contact structures by means of characteristic foliations is most 
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fruitful. Given a two-dimensional surface S embedded in M, the character- 
istic foliation of S, S%, is the [singular] one-dimensional foliation generated 
by the intersections of the tangent planes of T p S with the contact planes £ p in 
T p M. For all intents and purposes, may be thought of as a vector field on 
S generated by £ (by orienting the foliation). The singularities which arise 
on a characteristic foliation are generically saddles or spiral sources/sinks: 
pure centers cannot ever appear from a contact structure — the plane field 
twists too much for this. 




Figure 1. A contact structure is overtwisted if the induced 
characteristic foliation on some embedded disc possesses a 
limit cycle. 

The dynamical properties of S% are closely related to the topological classi- 
fication of contact structures. A contact structure £ is said to be overtwisted 
if there exists an embedded disc D C M such that D% possesses a limit cycle 
- a closed orbit along which nearby orbits accumulate [see Figure |lj. A 
contact structure is said to be tight if there are no such overtwisted discs 
anywhere in M. The contact structures for dz + xdy and dz + r 2 d9 are 
tight 0], while that of cos rdz + r sin r d9 is overtwisted {e.g., at the disc 
{r < 1, z = r 2 }). 

It is by no means apparent that the above definition is at all helpful: but 
in fact, the entire topological theory hangs on this dichotomy. Many major 
questions about contact structures are solved in the overtwisted category and 
unknown in the tight category (or, if known, then only recently and then 
by great skill and effort). For example, while overtwisted contact structures 
have been completely classified up to isotopy ||, the classification of tight 
structures appears [from what is presently known] to be delicate at best, 
intractable at worst |4], 
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4. Contact dynamics 



The connection between contact structures and Beltrami fields is, in the 
present context, quite straightforward. Given any fixed-point-free Beltrami 
field X on M, the Beltrami condition states that 

(4.1) MAX, •, •) = d(g(X, •)). 

^From this one can derive the crucial observation that the plane field or- 
thogonal to any nonvanishing Beltrami field is indeed a contact structure 
as follows. Witness the 1-form a := g(X, •) dual to X via g. The kernel of 
this 1-form represents the orthogonal plane field £ to X. This form a is a 
contact form on M since 



(4.2) a Ada := Xg(X, •) A[i(X,-,-), 

which, for A ^ 0, is nowhere vanishing, as one can easily check by evaluating 
on local orthogonal coordinate bases of the form (ei := Xj \\X\\ , 62, e^). 

A little more is in fact true: a Beltrami field annihilates the exterior deriv- 
ative of the associated contact form, since 

(4.3) (da)(X,-)=n(X,X,-)=0. 

Such vector fields are classical objects known as Reeb fields. The Reeb field 
of a contact form a is the unique vector field Z such that da(Z, •) = and 
a(Z) = 1. We have thus observed that any Beltrami field X (nonsingular 
with nonzero eigenvalue) is a Reeb field for a contact form, after a possible 
rescaling to force a(X) = 1. In jyj, a broader version of this result was 
demonstrated: namely, that the class of nonsingular Beltrami fields on M 
(up to scaling, for any Riemannian structure and volume form) is identical 
to the class of Reeb fields (up to scaling, for any contact form) . 

This theorem allows one to build "custom" solutions to the steady Euler 
equations of very high regularity. For example, |]l^] builds a single Beltrami 
field on a Riemannian M 3 which possesses closed orbits of all possible knot 
and link types simultaneously. This can be viewed as a rigorous manifesta- 
tion of the delightful results of Mofiatt on knotting in the Euler equations 



In the present context, we will use this simple correspondence between Bel- 
trami and Reeb fields to import technology from contact dynamics. Most 
specifically, we are interested in the utility of the tight/overtwisted di- 
chotomy in describing the dynamics of Beltrami fields. One extremely im- 



portant result in contact dynamics is the following theorem of Hofer [23]: Let 
£ be an overtwisted contact structure on M a compact 3-manifold without 
boundary. Then the Reeb field of any contact form associated to £ possesses 
a closed orbit, some multiple of which bounds a disc in M. 
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Loosely speaking, a closed orbit (limit cycle) in the characteristic foliation of 
a disc in M implies the existence of a closed Reeb orbit which bounds a disc 
in M. The proof relies on the delicate techniques of pseudo-holomorphic 
curves in symplectic manifolds. One indication of the implicit nature of the 
proof is that the location of the overtwisted disc has little to no correlation 
with the location of the implicated Reeb orbit. 

It follows from well-known properties of pseudo-holomorphic curves that the 
proof of Hofer's theorem remains valid for a three-manifold with invariant 



boundary (see [14] for details). Thus, for the solid torus, it follows that an 
overtwisted Reeb field must possess a contractible periodic orbit (since the 
fundamental group of the solid torus contains no elements of finite order 
except the identity). 



5. Definition of I 

^From Hofer's theorem, then, one way to force a contractible orbit in a 
Beltrami field is by finding an overtwisted disc in the orthogonal contact 
structure. This is far from trivial, since it requires searching for overtwisted 
discs among all possible embedded discs in M: not a computationally fea- 
sible task, even if the vector field were known analytically (which, in the 
context of an experimentally generated flow is not generally the case). 

The classification of contact structures has been successfully completed only 
on a selected class of three-manifolds. The classification of contact struc- 
tures on the solid torus is quite recent and subtle [|^, |24]| . In particular, it 
is known that on the solid torus there are tight contact structures which are 
"stably overtwisted" — taking some finite covering space of the solid torus 
and lifting the tight contact structure downstairs yields an overtwisted con- 
tact structure on the cover |14| (see Theorem |5.2| ). While this is a 
complication for contact topologists, it is a benefit to dynamicists. 

Lemma 5.1. Any Beltrami field transverse to a tight contact structure, 
some cover of which is overtwisted, must possess a contractible closed or- 
bit. 

Proof: Assume that £ is a tight contact structure, some finite cover £ of 
which is overtwisted. Then, given any Beltrami field X associated to £, 
lift this to a Beltrami field X on the overtwisted cover. Applying Hofer's 
theorem to the cover implies the existence of a contractible periodic orbit 
for X; however, since a covering space projection takes orbits to orbits, the 
closed orbit upstairs (along with the disc that it bounds) must project to a 
contractible closed orbit of the original Beltrami field X. □ 

Thus, our goal is to effectively determine the existence of an overtwisted 
or stably (with respect to coverings) overtwisted structure on a solid torus 
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given the least amount of information about a Beltrami field transverse to 
it. 

To do so, we recall a common index used in contact topology (see [Q, ||] for 
an introduction). Given a contact structure £ on a three-manifolds M, a 
simple closed curve (knot) is called transverse if its tangents are everywhere 
transverse to the contact planes. Assume that 7 is an oriented simple closed 
curve in M which bounds a compact oriented surface X in M. Then the self- 
Unking number of 7 with respect to £ and £ is defined as follows. Choose any 
vector field Z on a neighborhood of £ which has no fixed points and which is 
always tangent to £. That this is possible is a simple argument involving the 
classification of plane bundles. Then, flow 7 for a small amount of time under 
Z to obtain a "push-off" curve ■jz- The self-linking number of 7, sik^) is 
then defined as the intersection number of jz with X — i.e., the number of 
transverse intersections, counted algebraically using the orientations. This 
integer, which can be shown to be independent of the vector field Z chosen, 
is an invariant of transverse curves up to isotopy through transverse curves 
On S 3 , the self- linking number is also independent of the surface £ 
chosen so long as it bounds 7. 

The following recent result allows for an application of this index to Beltrami 
fields. 

Theorem 5.2 ((14]]). Assume a is a positive contact form on a solid torus 
V whose Reeb field is tangent to the boundary dV . Choose any transverse 
curve 7 on the boundary torus dV which bounds a meridional disc inV . If 
the self-linking number slk{^) of this meridian is not equal to —1, then the 
pullback of a under some finite cover is an overtwisted contact form. 

It known from the inequality of that if the initial contact structure is 
tight, the self linking number must satisfy slk < —\{D) = — 1. Thus, any 
self- linking number greater than —1 automatically implies an overtwisted 
structure (which is of course preserved under covers). The nontrivial re- 
sult of this theorem is that for a tight structure, a self-linking number less 
than —1 implies an overtwisted cover. The techniques used in the proof 
of this theorem are a combination of perturbing characteristic foliations, 
manipulating singularities of characteristic foliations, and using dynamical 
properties of the characteristic foliation on dV . 

Remark 5.3. It is necessary to distinguish between positive and negative 
contact structures. In the case of a negative contact form (one for which 
aAda < 0) on an invariant solid torus, the structure possesses an overtwisted 
cover if and only if the self-linking number of a transverse meridian is not 
equal to +1. It is an easy exercise to show that the sign of the contact form 
dual to a curl eigenfield is precisely the sign of the eigenvalue. 

^From these ingredients the following index may be defined: 
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Definition 5.4. Given a nonsingular Beltrami field X on an invariant Rie- 
mannian solid torus V, define the index I as follows. 

1. If the eigenvalue A of X with respect to the curl operator is zero, define 
I := 0. 

2. Otherwise, consider the characteristic foliation (dV)^ of the contact 
structure orthogonal to X on dV. If possible, choose 7 any meridional 
curve on dV transverse to dV^ and define 

(5.1) I := Sign(A) (s£k(j)) + 1. 

3. If no transverse curve 7 exists, define I := Sign(A). 

Theorem 5.5. Any C 2 or smoother nonvanishing Beltrami field on an in- 
variant Riemannian solid torus V having nonzero index I possesses a con- 
tractible closed orbit. 

Proof: Assume that A > as the negative case follows similarly. Since 
I 7^ 0, we are either in the case where there is a transverse meridian on 
the boundary of V with self- linking number not equal to —1, or there is 
no transverse meridian. In the case where the transverse meridian exists, 
some finite cover of the Beltrami field has a contractible closed orbit which 
is preserved by the covering projection. 




Figure 2. A Reeb component in a two-dimensional folia- 
tion is a foliation of an annulus by curves which limit onto 
the boundary components as illustrated [left, identify top and 
bottom]. On the boundary torus dV, if the characteristic fo- 
liation (dV)^ possesses a Reeb component [right], then there 
does not exist a transverse meridional curve. 

If there does not exist a closed transversal, then a basic result in foliation 
theory implies that either (1) the characteristic foliation (dV)^ is entirely by 
meridional curves; or (2) the characteristic foliation (dV)^ possesses a Reeb 
component, illustrated in Figure |2[ In the former case, the contact structure 
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is clearly overtwisted: any meridional disc in general position which spans 
one of these meridional curves has this boundary curve as a limit cycle in 
its characteristic foliation. The existence of a periodic orbit then follows as 
earlier. 

In the latter case, where a Reeb component exists in (dV)^, the following 
argument eliminates this possibility as a Beltrami field. Consider the Reeb 
field Z associated to the contact form a dual to X. Since Z preserves the 
contact structure £ and the boundary tours dV, it must likewise preserve 
the characteristic foliation (dV)^. Thus, if (dV)% contains a closed curve, 
then, since Z is everywhere transverse to (dV)e, the entire boundary torus 
is swept out by forward images of this curve under the flow of Z, and (dV)^ 
is a foliation by closed curves. Thus, a Reeb component (which always 
possesses both closed and open curves as in the illustration) cannot arise as 
the characteristic foliation on an invariant solid torus. □ 



6. Computation of I 

Given a contact structure £ on a 3-manifold M, the determination of whether 
it is a tight or overtwisted structure is a difficult question in general. By the 
Darboux Theorem (see, e.g., Q ^]), every contact structure in dimension 
three is locally equivalent to the kernel of dz+x dy, which is a tight structure 
[||. Thus, on the one hand, the property of being overtwisted is a decidedly 
global feature. However, the process of Lutz twisting allows one to 
change a tight structure into an overtwisted structure by means of a C° 
alteration on an arbitrarily small open set in M. Thus, given a Beltrami 
field X on M, determining whether the contact structure orthogonal to 
X is overtwisted is computationally intractable. Determining whether the 
universal cover is overtwisted is no less difficult. 

However, to compute the index I, one does not need information about 
the vector field on the entire three-dimensional regime, but rather on some 
(arbitrary) two-dimensional meridional disc. We outline a method for easily 
computing I from a C 1 approximation to X along a two-dimensional slice 
of the solid torus. 

Choose a meridional disc D with boundary curve 7. Orient 7 so that the 
contact form a := g(X, •) evaluates to a positive number on the tangents to 
7: i.e., 7 points in roughly the same direction as X. This orientation on 7 
induces in the usual way an orientation on the disc D. Denote by Xd the 
projection of the vector field X on to D by orthogonal projection onto the 
tangent planes (orthogonal with respect to the metric g). 

Proposition 6.1. If the vector field Xd is generic (possesses a finite num- 
ber of nondegenerate rest points), the index I of the Beltrami field X can be 
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computed by 

(6.1) I = Sign(A) 1 + Yl °(p)^d(X D ;p) 

\ p:X D (p)=0 J 

where for every rest point p of Xq, o~(p) is defined to be the sign (+/— ) of 
the dot product of X(p) with the positive normal vector to D at p, and the 
term Ind{Xo;p) denotes the standard Euler-F 'oincare index of the planar 
vector field Xjj at p. 



Proof: Assume that the characteristic foliation is generic in the above 
sense: there are a finite number of singular points p at which the contact 
structure £ is tangent to D, and the characteristic foliation about these 
points appears locally as a source/sink or a saddle. Following [Q, [J, there is 
a standard formula for computing the self-linking number of the transverse 
curve 7 = dD : 

(6.2) s£k( 1 )= J2 sign(T pJ D,e p )Ind( J D e ;p) 

Here, the sign of (T p D, £ p ) is +1 when the orientations on the contact plane 
£ and the orientations on the tangent plane to disc D at p agree. Otherwise 
the sign is —1. However, we only know the Beltrami field X and the disc D 
— not the characteristic foliation. Determining the indices of the rest points 
of the characteristic foliation is accomplished via the projected field Xd 
as follows. Since the contact structure for X is the plane field £ orthogonal 
to X, the characteristic foliation at every point q £ D is given by 

Dt{q) :=Z q nT q D=(X D (q)) ± , 

the line field orthogonal to Xjj at q. The proposition is proved by noting (1) 
fixed points of Xd occur exactly at fixed points of D^; (2) the Euler-Poincare 
index of at a fixed point p is unchanged by looking at the orthogonal 
vector field, as illustrated in Figure ||. □ 

Computationally, this is extremely simple as D can be chosen almost ar- 
bitrarily (one presumably chooses a D which is "nice" in coordinates) and 
information about X is required only on D itself. The local index calcu- 
lation is the most delicate portion of the computation: the location of the 
orthogonal point p is easy and the sign o~(p) merely measures whether X 
agrees with the oriented normal to D, which is trivial to determine. 

An example of a characteristic foliation and the resulting self-linking 
number is given in Figure ||. 
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Figure 3. Taking the orthogonal line field (grey) to the 
characteristic foliation (black) leaves the index invariant. 




Figure 4. An example of an sik computation given the 
characteristic foliation on an oriented spanning disc D. [left] 
D^; [right] the directions of the field X, positive to the 
left and negative to the right. The self-linking number is 
sik = — 3. Thus the Beltrami field X has index I = — 2. 

7. Conclusions 

The index I is an unusual object in that one inputs information about the 
vector field which is strictly iwo-dimensional, yet one obtains data about the 
dynamics which is fully i/iree-dimensional. In Figure ||, the only information 
about the Beltrami field known is that (1) it is orthogonal to the disc D 
at five points; (2) locally near those five points the projected field is either 
source/sink or saddle type; and (3) at those five points the field points out or 
in as illustrated. From this data, it is inevitable that somewhere in the flow 
there exists a contractible closed flowline. This is a corollary of the contact- 
topological methods used — the moral of the story is that the Beltrami 
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condition hides within it certain constraints on the dynamics which couple 
the dynamics of the vector field to the topology of the orthogonal plane field. 

A deficiency of our theory is that it is not sharp. It is certainly possible 
for a Beltrami field to have the value 1 = and yet still have a contractible 
periodic orbit: indeed, any tight contact structure which has been Lutz 
twisted along a contractible closed curve necessarily has trivial index as 
well as a contractible orbit. It appears certain (due to this mechanism of 
arbitrarily small Lutz twists) that no completely sharp computable index 
can be defined. What I does, however, is detect if there are contractible 
orbits forced by the presence of non-localized overtwisted discs, and in this 
regime it is efficacious. 

It would be interesting to find a sharp lower bound on the number of periodic 
orbits present in the case of I 7^ (cf. the recent body of theory surrounding 



the contact homology of Eliashberg, Givental, and Hofer [10Q. 
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